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Abstract 

We consider a class of higher order corrections in the form of Eu- 
ler densities of arbitrary rank n to the standard gravity action in D 
dimensions. We have previously shown that this class of corrections 
allows for domain wall solutions despite the presence of higher powers 
of the curvature. In the present paper we explicitly solve the lin- 
earized equation of motion for gravity fluctuations around the domain 
wall background and show that there always exist one massless state 
(graviton) propagating on the wall and a continuous tower of massive 
states propagating in the bulk. 
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It may very well be that the theory of gravity as we know it today is 
only an effective theory and the usual Einstein-Hilbert action should be sup- 
plemented with corrections involving higher powers of the curvature tensor. 
This point of view is supported for example by string theory or the presence 
of the conformal anomalies in all quantum field theories coupled to gravity. 

In this paper we are interested in the corrections of special type - Euler 
densities of arbitrary order n (n being a power of the curvature tensor) in 
arbitrary space-time dimension D. It was shown in p[ that in the presence 
of arbitrary number of Euler densities in the lagrangian, there always exist 
domain wall solutions. The order n should be less or equal to D/2 where D 
is the dimension of space-time {D = 2n although formally total derivative 
gives part of the conformal anomaly). Euler densities appear for example in 
the a' expansion of the string theory effective action 0. The Euler density 
of the order n = 2 (equal in this case to the Gauss-Bonnet combination) has 
been discussed in the presence of branes in many papers, some of which are 

i-ii- 

In the present paper we analyze the linearized equations of motion for 
the fluctuations of the metric around the domain wall solutions in the theory 
with arbitrary number of Euler densities. It is a generalizations of the original 
idea [jl5| and the analysis performed for Gauss-Bonnet n = 2 case 0] (the 
Gauss-Bonnet term for intersecting domain walls has been recently discussed 
in [|14|)- The fluctuations are assumed to be graviton-like i.e. transverse and 
traceless. 

It turns out that in spite of the presence of higher powers of curvature the 



picture is very similar to the usual Randall-Sundrum scenario . There ex- 
ists one normalizable massless bound state and a continuous tower of massive 
states with small amplitude on the wall. The strategy adopted to calculate 
the equations of motions is based on the explicit formulae for the Euler den- 
sities derived in our previous paper [0. 

The metric of the domain wall is conformally flat (and this fact was exten- 
sively used in 0) but this is no longer the case when we add a fluctuation to 
the metric. Therefore calculating the equations of motion for the fluctuation 
requires calculating a first order correction (linear in the Weyl tensor) to the 
formulae in 0]. When adding all the contributions, the resulting equation 
of motion turns out to be almost identical to the lowest order case, the only 
difference being in the actual value of the coefficients of the equations. It 
is extremely important to notice that the graviton (normalizable massless 
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mode confined to the wall) exists independently of the number of dimensions 
and presence of the Euler densities of higher order. It is quite amazing that 
the picture in presence of Euler densities is almost identical to the lowest 



order scenario 15 



The Euler densities in D dimensions are defined (in the form notation) 



as 



= ^ rreaia^-ani?""'' A ■ ■ ■ A A e''^"+' A ■ ■ ■ A 6^° . (1) 

{D-2n)\ ' ' ^ ^ ^ 

(For D = 2n they are topological invariants and formally total derivatives 
but a careful regularization shows that one cannot discard them neither in the 
action nor in the equations of motion since they correspond to the conformal 
anomaly) . 

We will consider models in which D-dimensional gravitational interactions 
are described by the sum of such Euler densities and in which there is a 
{D — 2)-brane (a domain wall). The action is the sum of the bulk and brane 
contributions: 

S = "S'bulk + 'S'brane 



5 

d'^xy^ ^ , (2) 

n=0 

^brane = / d^'^X^Q (-A + X^R + . . ) . 



The metric on the brane is given by g^i,{xP) = g^u{xP,y = 0) where y = x^ 
and /i, z/, . . . = 1, . . . , D — 1 while M, N, . . . = 1, . . . , D. In the brane action 
we write explicitly only the most important term - the brane cosmological 
constant. The first two terms of the bulk action are known from conventional 
gravity. The one with n = corresponds to the cosmological constant: 
/(o) = 1, = -A. The one with n = 1 is the usual Hilbert-Einstein term, 
I^^^ = R with the coefficient ki = (2/t^)~^. The maximal number of the 
higher order terms is rimax < [-D/2] as discussed in [0]. 

We assume in the following that Ai = 0; Ai 7^ would give additional 
contribution Aim^/im.(0)/2 to the right hand side of (^) discussed later so it 
would affect only the massive modes. 

Let us start with the bulk equations of motion. They can be obtained 
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from the variation of the vielbein in the bulk action (for n < [D — l)/2) 

E ""p^ ~^jj^ ga,a....a^_,ai?"^"^ A- ■ ■ Ai?"^-^"^" Ae"^"+^ A ■ ■ -Ae""'^ = 0. (3) 
We can write the curvature two-form as 

^ah ^ (jab ^ _1_ A - 6^ A K"") (4) 

where C"'^ is a two-form composed of the Weyl tensor Cj\/7v_r5 while i^'" is 
a one-form which will play important role in our calculations and which is 
defined (for invertible vielbeins e1[) as K'^ = KuN^^'^dx^ = K"-i,^^ with 

Kmn = Rain — —^^QmnR ■ (5) 

The main purpose of this paper is to analyse localization of the effective 
brane gravity in higher order theories proposed in our previous paper [0]. 
In order to do this we will look for solutions of the linearized equations of 
motion (|]) for the fluctuations of the metric: 

dMN = dMN + ^^MN (6) 

for which the background metric Qmn is of the domain wall type that was 
proven in @ to satisfy the equations of motion (|]). The line element of this 
domain wall background is equal 

ds^ = e-^f^y^^^dx^'dx" + dy^ (7) 

with the warp factor function given by f{y) = <j\y\- The fluctuation Hmn is 
assumed to propagate along the brane and be transverse and traceless: 

hoM = 0, v^'h^u = k^Ku = . (8) 

We decompose the general fluctuation into modes with definite mass (from 
the [D — 1) -dimensional point of view) 

h,A^'',y)=i:e''-^'h^M (9) 
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where k is a. {D — 1) -dimensional momentum satisfying k"^ = —m? and the 
sum is for discrete modes while the integral is for continuous modes (in the 
above formula and in the next one the factor exp(i/Co-a;'^) represents the sum 
of independent real fluctuations with a given mass). Let us now concentrate 
on one mode hm^^{y)- Calculating the curvature tensor for the metric 

gDD{x'',y)=g%{x%y), (lo) 

the equation of motion will be expanded to the first order in e. For the 
background metric (|^) the Weyl tensor vanishes therefore C"* is already of 
order e: 

= e{CiT'' + C(e') . (11) 

Thus, only terms up to the first order in (Ci)"'' should be kept in (|^). There 
is a second contribution to the linearized equations of motion for h^^^ coming 
from the correction to K"^, 

= (Ko)" + e(^i)" + Oie^) , (12) 

and in the product (^ also for this contribution only the first power of (-ft'i)" 
should be kept. 

Using the above e expansion we can write the equation of motion (^) up 
to terms linear in e in the form: 

° = ^ ? ''" (D-mD-l-2ny. ((^«"''« + '(^''"')« + ''^'"'>") 

where we introduced the following tensors (note that Hq^'* have different 
normalization than the corresponding tensors i/*^") of ref. [0) 

= K/^m/S ■ ■ ■ ± perm.) (i^^o)K^(i^o)t " " " (^o)i^: , (M) 
= n(«;4^■■^^:±perm.)(KOS(iro)^^..(i^o)t, (15) 

/ („)xiv ^ n{D - 2) 

■ ■ ■ ■ ± perm.) (Ci)KH^o)t ■ ■ ■ (Ko)^: ■ (16) 
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In it was shown that the lowest order (e*^) term vanishes for the domain 
wall metric (0) with the warp factor o determined in terms of the coupling 
constants Kn (this domain wall is flat when some relation among k„ is satis- 
fied). 

Let us start the discussion of corrections of order e from . From (|15|) 
one can get the recurrence relations for h'i^ (valid not only for the domain 
wall but for all background metrics): 



)p 



(17) 

To proceed further let us specialize to the transverse and traceless fluctua- 
tions of the metric (^) in the domain wall background (|^). In such a case the 
tensor Ki satisfies the following conditions 



(^i)i: = 0, (K,)^ = {K.Yd 



Tensor Kq has been calculated in P 



(18) 



(^( 







)d 



,. ^-2 (df 
' 2 [dy^ 

D-2 fdfV 



+ 



dy J dy 



while Hq""^ differ only by normalization from H^^'> defined in 



(n)w 



D 



{D - 1 - n)! 



2n-2 



dy, 



2n fd^f 



D-1 \dy^ 



(^. 



)d 



iD-l-n)\ [d^^ 



2n 



(19) 



(20) 



Substituting eqs. (|r^) and ( PUj ) into formula ( P!7[ ) and using the conditions 
(ITSD we get the following expression for h[^^ 



(21) 
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with all other components vanishing. It explicitly reads 

'D-2\"-^ n{D -2)1 




2n-4 




1 -n)! 

2 



2(n-l) /d'f 
D-2 \dy^ 



(22) 



Let us now turn to G\ \ For a general background metric g%]\f the 
expression obtained by taking into account all permutations in eq. (|I6|) is 

(k) 

very complicated. The result contains different combinations of all Hq with 
k < n and it is not possible to write it as a simple recurrence analogous to eq. 
(|T^ valid for h[^\ Therefore, we will present explicit formulae for G^^^ only 
for the domain wall background (J^). In such a case, using the symmetry 
properties of the Weyl tensor C and the fact that (Kq)'^ and [h'q^Y^^ 
proportional to 5^, we get 



are 



YD n\{D 



E 

k=l 



D~l-n 



[n-k)\{D-iy 



0)p 



k-l 



{n-k)^ 



(23) 



Substituting the explicit formulae for Kq and H^'^^ (|T^ED|) it can be rewritten 
in the following form 



(n 



{D-1 



n] 



D 




2n-4 



(24) 

We see that the tensors ifj"^ and G^i'' are proportional to Ki and Ci, 
respectively, wich can be found using the decomposition of the curvature 
tensor (given by eq. 




2 2f 

-m e 



h^l{y). (25) 
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D-3 
D-2 



1 92 



d 



2dy'^ 



2 \dy ) dy 

2 / 02 . 



9y\ 



2{D-3) 



2 2f 

m e ■' 



hm'Jy). (26) 



In both these tensors the operators acting on h^'^^^y) do not depend on the 
Lorentz indices n, v so we will drop these indeces from now on. 

Now we are ready to get the equation of motion for the fluctuation hm{y)- 
We rewrite equation (0) using the explicite form of h["'\ G^i \ Ki and Ci 
eqs. (^ , p4| -p^) ■ The part of this equation linear in e gives: 







E 



Kr. 



2n{D - 3)! 



2n)\ 



2 ^ dy^ 



D 

+ — 



+ {D - 3) {fT - i^n - 1) (/') 



dy 



d_ 

dy 



if) - '-m\y UT-' 



n 



D-3 



m e ^ [f) (/ ) 



hm{y) 



(27) 



At this point we can use the explicit form of the warp factor f{y) = a\y\. 
One should be careful when performing this substitution because the first 
derivative of f{y) is not continuous at ?/ = and the second derivative of 
f{y) is proportional to the Dirac delta at y = 0. Similar care is needed when 
calculating terms containing derivatives of the fluctuation because the 
solutions to the above equation are also functions of \y\. Regularizing the 
Dirac delta function one can find the following equalities^: 



(/ 



>\2k 



ifTU") 

:hm{y) 



2k + I 

^2ku" 



-a 



2fc+l 



dy' 
d 



a 



2k + 1 



a''S{y)h'^{0' 



ifT-\n^hUy) 



2k 



(28) 
(29) 
(30) 



^ The necessary denominators of 2fc + 1 have not been taken into account for example 
by the authors of ref. Q (who discuss the case D — 5) leading to wrong equations of 
motion. 



8 



where 



(31) 



Making these substitutions in eq. ( P7D we find the following bulk (y 7^ 0) 
equation of motion for hm '■ 







n{D-3)\ 



2n-2 



' {D-1- 2n)\ 
[-hi + {D- 5)ah'^ + 2{D - 3)a^hm - m~e 



(32) 



It is rather amazing that the expression in the square parenthesis does not 
depend on n. Therefore for arbitrary k„ and Umax the bulk equation of 
motion for hm reduces to 



= -h'l + {D- 5)ah'^ + 2{D - 3)a^hr 
Its solution for m = is equal to 



m' 



(33) 



AD~'i)^\y\^ (34) 
while for m 7^ it can be written using the Bessel functions of order 



Arr,.J D-1 



m 



a 



m 



a 



(35) 



For each m only one of the above combinations of the solutions of the bulk 
equation ( p^ ) is the solution of the full equation (^). In order to identify 
this combination we have to take into account the part of the equation of 
motion (p^) proportional to S{y). It reads: 



2n) 
2h^{0) 



■a 



2n-l 



2(n-r 



m 



/im(0) 



a ■ • (D-3)(2n-3) or2 

This equation is equivalent to the following boundary condition at y = 0: 



(36) 



hm{0) 



-2a 1 



m 



(D-l-2n)! (2n-3) 



CJ^D-3 Y.n^^n{-iy 



(D-l-2n)! 



cr 



2n 



(37) 
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For m 7^ it is quite complicated and depends on k„ but in principle can be 
used to fix the ratio of the coefficients Am/Bm for the massive modes. 

The situation is much simpler for m = because in this case the above 
boundary condition does not depend on k„ and simplifies to 

Applying this condition to (^) we find that for arbitrary space-time dimen- 
sion, D, and for arbitrary strength of the higher order interactions (given by 
the coefficients k„) the massless solution of the equations of motion always 
exists and is given by 

ho{y)=exp{-2a\y\). (39) 

We would like to interpret these solutions as massless, normalizable 4- 
dimensional gravitons. It turns out to be possible, but some care is needed. 
Let us start with the normalizability issue. To check whether the above 
solutions (massless and massive) are normalizable we have to choose an ap- 
propriate integration measure. Such measure can be determined by the re- 
quirement that the kinetic operator for the gravitons should be self-adjoint 
(which is not fulfilled in the case of eq. (|32D). To find the proper operator it 
is helpful to change the variable y and to rescale fluctuation hm- 

hm{y) = {I + a\z\)^ hm{z) (40) 

where 



az = sgn(y) (e'^l^l - 1 j . (41) 
The equation of motion reads then 

d\ , , D{D - 2)a\ , , 2, . N /.^N 

- -T^hm{z) + — — —r-:rhm{z) = Vfl hm{z) (42) 

dz^ 4(l + cr|z|)^ 

and is explicitly self-adjoint with a fiat measure (which we can take equal to 
!)• 

The properly normalized massless solution is given by 



ho{z) = /-^-^ {1 + a\z\)-^ (43) 



D-2 
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and is normalizable (for D > 3) while the massive modes are given by the 
formula 



hm{z) = \/—+m\z\- (44) 
AmJ 1- rn\z\ + Bm — Y d-i h m\z\ 

Massive modes are asymptotically (for large \z\) plane waves and therefore 
for infinite range of z not normalizable. 

We showed that the massless mode is always normalizable but this is not 
enough to interpret it as a graviton. The reason is the following. The bulk 
equation of motion ( 127|) contains an overall factor depending on k„ and on 
the warp factor a: 

i^"^"^ (D-l-2n)!^ • ^^^^ 

Although the value of this factor is not important for the solution, we have 
to remember that the sign of this factor is related to the sign of the kinetic 
energy term for the fluctuations in the effective lagrangian. The wrong sign 
of the kinetic energy indicates instability of the assumed background. Thus, 
our domain wall solution can be stable only if the parameters such 
that the sum in eq. (^) is positive. Using rescaled couplings (introduced in 

i): 

= ^-i)n-i^ (D-iy. 

yn V ; "(D-l-2n)! ^ ^ 

the necessary condition for the stability of the domain wall solutions can be 
therefore written in the form 

Y: np^a'"^-' > . (47) 

n=l 

It is interesting to compare the above condition with the bulk and bound- 
ary equations which must be satisfied by the warp factor, cr, of the back- 
ground domain wall metric (|^: 



J2p^a'^ = -A, (48) 



n=l 



E ^^Pna^"-^ = ^-^\ . (49) 
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In our previous paper 0] we have considered possibility of the domain wall 
solutions without the bulk or/and the brane cosmological constants. Such 
solutions are acceptable only if the new condition (^) is satisfied and this 
must be checked for any specific model described by a set of the coefficients 
It is not easy to discuss the consequences of (|47|) in general but one can 
make the following observation. 

Let us assume that we insist on A = - the "self-supporting" brane 
solution without any matter on the brane (the bulk cosmological constant 
may be different from 0). Let us define a polynomial 

n^) = E ■ (50) 

n=l ^'^ ^ 

The condition for vanishing A reads -P(cr) = while the condition for the 
correct sign of the kinetic terms is P'{(j) > 0. From the positivity of the 
gravitational constant ki = {2k^)~^ it follows that -P'(O) > 0. This means of 
course that the trivial domain wall with a = (i.e. the Minkowski space) 
has the correct sign of the graviton kinetic term. This means also that the 
first nontrivial domain wall (the one with the smallest positive a) with A = 
is unstable. But not all solutions with A = must be unstable. The number 
of different solutions with A = and o" 7^ is equal to the number of positive 
zeros of P{x) which is at most {rimax ~ !)• If P{x) has only first order 
zeros than "every second" solution with A = can be stable because such 
function changes is derivative when moving from one zero point to the next 
one. Thus, the nontrivial domain wall with vanishing brane cosmological 
constant is possible if rimax ^ 3 and this can be satisfied for space-time 
dimension D > 6. 

Similar analysis can be performed for domain walls with vanishing bulk 
cosmological constant, A = - the only difference is to use the polynomial 
J2n=T Pnx'^'^ instead of that defined in eq. (^. 

In the case of both cosmological constants vanishing it is not difficult 
to see that for rimax > 3 there exists a range of values of pi, . . . ,Pnmax ^^r 
which all three conditions (47-^91) can be simultaneously satisfied with A = 



A = and a ^ 0. Thus, stable solutions with vanishing brane and/or bulk 
cosmological constant are possible if the space-time is at least 6-dimensional. 

In conclusion, we considered a class of models with higher order grav- 
ity corrections in the form of the Euler densities with arbitrary power n of 
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the curvature tensor in arbitrary space-time dimension D. The fluctuations 
around the domain wall type solutions (found in 0) were shown to have sim- 
ilar spectrum as in the lowest order case (n = 1) - the bulk equation of mo- 
tion rather miraculously turned out to depend only on D and not on n. The 
boundary condition at the wall for massive modes has some n-dependence. 
There exists one normalizable massless mode and a continuum of massive 
modes (without the energy gap). The solutions for all D are almost the 
same as in the original Randall-Sundrum model with the Hilbert-Einstein 
action |jl5| (apart from some numerical factors) and the discussion about 
the applicability of the Newton's law and the effective number of dimensions 
can be carried over from ([|15|) to the general case discussed in this paper 
virtually unchanged. 



This work was partially supported by the Polish KBN grants 5 P03B 150 
20, 2 P03B 052 16 and by the European Commission RTN grant HPRN- 
CT-2000-00152. 



References 

[1] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[2] K.A. Meissner, M. Olechowski, Phys. Rev. Lett. 86, 3708 (2001). 

[3] R.R. Metsaev and A.A. Tseytlin, Nucl. Phys. B293, 385 (1987); K.A. 
Meissner, Phys. Lett. B 392, 298 (1997). 

[4] J. Kim, B Kyae and H. Lee, Nucl. Phys. B582,296 (2000); Erratum 
Nucl. Phys. B591, 587 (2000). 

[5] I. Low and A. Zee, Nucl. Phys. B585, 395 (2000). 

[6] Z. Kakushadze, Nucl. Phys. B589, 75 (2000); O. Corradini and Z. 
Kakushadze, Phys. Lett. B494 302 (2000). 

[7] M. Giovannini, Phys. Rev. D63, 085005 (2001); Phys. Rev. D63, 
064011 (2001). 

[8] N.E. Mavromatos and J. Rizos, Phys. Rev. D62, 124004 (2000). 



13 



[9] S. Nojiri and S.D. Odintsov, Phys. Lett. B484, 119 (2000); JHEP 07, 
049 (2000); |hep-th/0105lT7 . 

[10] I. P. Neupane, JHEP 09, 040 (2000); Y. M. Cho, I. P. Neupane and 
P. S. Wesson, |hep-th/01042"27| . 



[11] K. Kashima, ^p-th/001028e . 

[12] J. Kim, B. Kyae and H. Lee, |hep-th/0104T5q . 

[13] L. Anchordoqui, J. Edelstein, C. Nunez, S. P. Bergliaffa, M. Schvel- 
linger, M. Trobo and F. Zyserman, |hep-th/0106127 . 

[14] L P. Neupane, |hep-th/0106T00| ; 

[15] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 



14 



